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We show that a simple arrangement of external fields, consisting of a static component and 
an orthogonal rotating component, can continuously decouple a two-qubit entangled state from 
uncorrelated dephasing, bit flipping, and dissipation at finite temperature. We consider a situation 
where an entangled state shared between two non-interacting qubits is initially prepared and left 
evolve under the environmental perturbations and the protection of external fields. To illustrate the 
protection of the entanglement, we solve numerically a master equation in the Born approximation, 
considering independent boson fields at the same temperature coupled to the different error agents 
of each qubit. 
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I. INTRODUCTION 

Entanglement of qubits is a fundamental resource of 
quantum information processing A pure entangled 
state, initially prepared as a superposition of factorized 
states, inevitably decoheres under environmental noise 
and could even disentangle completely within a fi- 
nite time interval Q . Consequently, any potentially suc- 
cessful design of a device intended to perform quantum 
computation must anticipate ways to preserve qubit en- 
tanglement. 

There are very sophisticated methods developed to 
protect quantum states, like error-correcting codes Q 
and strategies based on decoherence-free subspaces and 
subsystems Q, where logical qubits do not necessar- 
ily coincide with physical qubits. This usually requires 
that the quantum information corresponding to N logi- 
cal qubits be stored in more than A'^ physical qubits. It is 
also possible to protect the quantum information stored 
directly in physical qubits by using the method of dy- 
namical decoupling . We expect that a fully functional 
quantum computer will use all these protecting tools in 
a coordinated way for maximum efficiency and fidelity. 
Here, due to its relative simplicity, we focus on the con- 
tinuous version of dynamical decoupling 0, Q and show 
that it can protect a two-qubit entangled state at finite 
temperature from uncorrelated dephasing, bit flipping, 
and dissipation. 

The paper is organized as follows. Section II presents 
the total Hamiltonian describing the qubit system and 
the rest of the universe. We also introduce the general 
method of dynamical decoupling employed to derive the 
simple control fields. The master equation for the evo- 
lution of the two-qubit density matrix is derived in Sec. 
III. The description of the environmental errors is given 
in Sec. IV, where we demonstrate the efficacy of the 
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continuous dynamical decoupling in protecting entangle- 
ment. Finally, we conclude in Sec. V. 

II. THE TOTAL HAMILTONIAN 

The Hamiltonian for the two-qubit system and the en- 
vironment is written as 

H{t)=Ho{t)+HE + H,nt, (1) 

where HQ{t) is the qubit-system Hamiltonian, including 
terms describing the action of the control fields, He is the 
environmental Hamiltonian, and iJjnt is the term repre- 
senting the interaction between the qubits and their sur- 
roundings. The general prescription for dynamical de- 
coupling ^, ^ consists of finding a control Hamiltonian, 
acting only on the qubit system, such that its correspond- 
ing unitary evolution operator Uc{t) is periodic and 

/ " Ul{t)H,^tUc{t)dt = 0, (2) 
Jo 

where tc is the period of Uc{t). In the following we ex- 
plain a possible strategy to find a combination of static 
and simple oscillating external fields leading to a periodic 
unitary operator Uc{t) satisfying Eq. ([2]). 

The interaction Hamiltonian iJjnt describing two 
qubits coupled to their surroundings is written as 

^int = Bi • (Ti + B2 • <T2, (3) 

where Bfe = J2m=i Bk,m^rn, for /c = 1,2, with xi = x, 
X2 = y, X3 = z, and Bk,m, for k — 1,2 and m = 1, 2, 3, 
are operators that act on the environmental Hilbert 
space. Here, for fc = 1,2, CTk = xak.x + ycrk,y + zo-fc,z, 
and Uk.x, '^k,y, and ak,z are the Pauli matrices acting 
on qubit fc. The unitary operator exp(—2inxTTtai^x /tc) 
is periodic with period tc for any integer Ux 7^ 0. If we 
identify this operator with Uc{t) in Eq. ^ and take iJmt 
as given by Eq. ([3]) , the integration does not give zero in 
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general, but eliminates all terms proportional to ai^y and 
ai^z- If, to replace Uc{t) in Eq. we consider the op- 
erator eyip{—2inx'^t(7i^x/tc) exp(— 2m27rtcri.2/fc) instead, 
where and are non-zero integers, then the integra- 
tion eliminates all terms proportional to cri if ^ Hz, 
there remaining only the term <cB2 •cr2. Hence, to satisfy 
Eq. (21), we choose 



Ucit)^U2it)Ui{t)^Ui{t)U2{t), 

since <Ti and cr2 commute, where 



Uk{t) = exp 



exp 



. 2n,7rt 



tr 



-crk,z 



(4) 



(5) 



for fc = 1,2. Hence, if the control Hamiltonian is writ- 
ten as Hc{t) = n{t) ■ {<Ti +(T2), then Eqs. (g]) and ^ 
imply the very simple external-field configuration of our 
previous work on singie-qubit operations, Ref. [s]: 



n{t) = 



n,uj z cos 



{n^ujt) - y sin [n^ut)] , (6) 



where uj — 27r/ic- 

To study the protection of entanglement against envi- 
ronmental sources of noise, we focus on a situation where 
the two qubits are prepared in an entangled pure state 
at time i = 0. We further assume that the qubits do not 
interact with each other and, if they could be isolated 
from the rest of the universe, their non-local state would 
remain unchanged. Thus, let us take Hq = ujo{'^i,z + o'2,z) 
as the unperturbed two-qubit Hamiltonian, using units 
for which h = 1. Let r be the time interval during which 
we intend to preserve the entanglement. We then choose 
tc = t/N, where N is an integer, so that C/fe(r) = /, for 
k = l,2. 

Equation ^ gives the Hamiltonian for the evolution of 
the ket representing the joint state of the qubits 

and their environments. For the sake of simplicity, we 
assume there is a global, static control field along the z 
axis chosen to cancel Hg exactly. Hence, the remaining 
terms of Hq (t) represent the action of additional control 
fields. We identify the evolution dictated by Ho{t) with 
the action of the unitary operator Uc{t) of Eqs. ^ and 

Here we show that the very simple field configuration of 
Eq. ([6]) can also prevent a two-qubit entangled state from 
disentangling due to uncorrelated dephasing, bit flipping, 
and dissipation at finite temperature. We emphasize that 
Eq. ([ni is a simple combination of a static field along the 
X axis and a rotating field in the yz plane. Moreover, ad- 
dressing each qubit independently is not necessary; the 
field is supposed to be spatially uniform in the neighbor- 
hood surrounding both qubits. 



III. THE MASTER EQUATION 

In the interaction picture, the Hamiltonian is given by 

2 3 

= uUt)Bk,^UE{t)Ul{t)<Jk^Mt), (7) 

k—1 m—1 



where ak,i = crk,x, (rk,2 = crk,y: o-k,3 = <^k,z: Usit) = 
exp{—iHEt), and we have used Eq. ([3]). The quantities 
U^{t)(TkjnUc{t), for A; = 1,2 and m = 1,2,3, are rota- 
tions of (7k,rm whose matrix elements, Rm.nit), are real 
functions of time: 



Ul{t)(7k,mUc{t) = Rm,n{t)(7k,n- 



(8) 



If we define the operators Ek,m{t) — U\,{t)Bk^mUE{t), 
for /c = 1,2 and m = 1,2,3, and use Eqs. ^ and ®, 
then the interaction Hamiltonian becomes 



Hl{l:) — ^ ^ Rm,n{t)Ek^m{t)(yk,n- 



(9) 



k—1 m—1 n—1 



In the interaction picture, the Redfield master equa- 
tion describing the temporal evolution of the two-qubit 
reduced density matrix, pi{t), is written as (lo| : 



dpijt) 
dt 



dt'TiE { [Hi {t) , [Hj {t') ,pEPim},m 



where we have assumed the noise is low enough that 
the Born approximation is valid. We also notice that 
Eq. (|10p is not a Markovian master equation, since 
the dynamical-decoupling process occurs in a time scale 
shorter than the environmental correlation time; this is 
the reason we keep t as the upper limit of the integral 
on the right-hand side of Eq. ^ (cf. p. 132 of Ref. 
[To|). Here, p_b is the initial environmental density ma- 
trix, pe = exjp{—(3HE)/Z , where Z is the partition func- 
tion given hy Z = Tte [exp(— /Ji/^;)], /? = l/fc^r, kB 
is the Boltzmann constant, and T is the absolute tem- 
perature, assumed to be the same in the surroundings 
of both qubits. By substituting Eq. ([9]) into Eq. (fTO]) 
we encounter the quantities Tte [Ek,mit)pEEk\m'it')], 
for k,k' = 1,2 and m,m' — 1,2,3. To illustrate our 
methodology in a simple manner, we suppose that the 
reservoir operators at the position of one qubit are un- 
correlated with the reservoir operators at the position 
of the other. Moreover, we also assume the qubits and 
their respective environments are identical. Thus, we 
can write Tte [Ek,„i{t)pEEk' ,m'it')] = 6k,k'Cm,m'it,t'), 
for k,k' = 1,2 and m,m' — 1,2,3, where Cm,m'{t,t') is 
the correlation function between components m and m' 
of environmental field operators calculated at the same 
qubit position. We, thus, define the quantities 
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E 

m,n— 1 



Rm.p{t) / dt' Rn,q{t')Cm,n{t,t'), (11) 




for p,q = 1,2,3. Hence, the master equation now be- 
comes 

dpiit) ' ' 



dt 



k=l p,q=l 
2 3 

+ E E Wk,gPi{t),<yk,p] , (12) 

k=l p,q=l 
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where we have assumed the environmental fields are Her- 
mitian. 



CONTINUOUS DYNAMICAL 
DECOUPLING 



IV. 



To solve Eq. ^ we need -Dp,g(i), Eq. ([IT]), and, 
therefore, we must calculate the correlation functions 
Cm,n{t,t'), for m,n — 1,2,3. Here, as stated above, we 
consider independent dephasing, bit flipping, and dissi- 
pation. Associated with these errors, we introduce six 
independent boson fields, three at each qubit position, 
all of them at the same finite temperature T. Accord- 
ingly, the terms appearing in Eq. ([3|) can be written as 

Bfc • (Tfe = ak,z ^ [si.Aflfe.A + giWk,\ 

A 



93,XCk,\ 



+ {<Jk,x - i<Jk^y) 9l\cl,\' 



(13) 



for k = 1,2, where 51, a, 52, a, and g^^x are complex 
coupling constants that do not depend on the qubit- 
position index k, reflecting the fact that the qubits are 
surrounded by identical environments, ak^x, b^^x, and 
Ck,x are, respectively, the annihilation operators for 
mode A of the boson field associated with dephasing, 
bit flipping, and dissipation, with respective creation 
operators a\ x^ ^1 A' ^^'^ '^k x- field operators depend 
on fc, since, although identical, the qubit environments 
are uncorrelated. The only non-zero commutators of 
these operators are: [ak,x,al^x\ = 1' \-^k,x,bl^^] = 1, and 

[cfc,A,c^A] = ^' k = 1,2, and all A. Therefore, we 
take the environmental Hamiltonian as given by He = 

ELi Ea [^hxal xak^x + i^2.xbl xbk,x + i^3.xcl xCk,x\ , 
where wi^x, W2,a, and tos^x are the Ath-mode frequencies 
of the fields associated with dephasing, bit flipping, 
and dissipation, respectively. Using Eq. we 
can calculate the correlation functions Cm,m' [t, t') — 

Tr£ [El^r,^{t)pEEl,m'{t')] = Trfi [E2^ra{t)pEE2,m'{t')], for 

m,m' = 1,2,3. The only non-zero correlations are: 





= /C2(^-^')^ 




t')] 






+ IC3{t-t') 


+ 2Re [C3it 




(14) 




= iK.3{t-t') 


-2Im [C3{t- 




(15) 


C2At,t') 


= ~ilC3{t-t' 


) + 2Im [C3ii 


~t')]. 


(16) 


C2,2{t,t') 


= lC3it-t')^ 


-2Re [£3{t~ 




(17) 




= JCiit-t')^ 


- 2Re [£1 {t ~ 


0], 


(18) 



where the complex functions JCmit) and Cm(t), for m = 
1, 2, 3, are given by ICm{t) = J2x \9rn,x\'^ exp{iu!rn,xt) and 

Cm{t) = J2x \9m.x\'^ exp{iu!rn,xt)/ [exp{f3uJm,x) " 1]- 




FIG. 1: A{t) for the evolution of different initial entangled 
states. The concurrence as a function of time is given by 
max{0, A(t)}. (a) A{t) for initial states with A(0) « 0.70711, 
when the control fields are turned off (indicated with A'^ = 0) 
and on (A'^ = 1,3,8). The solid and dotted lines correspond 
to the initial states of Eq. ([T9|l . with 6 = n/8 and 6 = Stt/S, 
respectively. The double-dot-dashed line corresponds to the 
initial states of Eq. ([20}, with S = 7r/8 and 6^ = 371/8 giving 
the same A{t). (b) A{t) for initial states with A(0) = 1 (BeU 
states). The initial states of Eq. (|19p evolve resulting in the 
same A(i) for 6 = ±7r/4. Also, the initial states of Eq. pO|l 
evolve resulting in the same A{t) for 9 = ±7r/4. Here we 
show the results for Eq. (|19|l . since they differ from those 
of Eq. (|20|l only in the fourth decimal place. The solid line 
represents the result for the control field turned off, while the 
dotted, double-dot-dashed, and dashed lines show the results 
for control fields with A'^ = 1,3,8, respectively. 



In the limit in which the number of environmental 
normal modes per unit frequency becomes infinite, we 
define spectral densities for m = 1,2,3 as Jm{(^) = 
J2x\9m,x\'^ 6{u! - ujjn,x), with u e [0,-|-oo) and inter- 
pret the summations in ICm{t) and Cm{t) as integrals 
over Lu: ICm{t) — Jm(a;) exp(ia;t) and Cm{t) — 

dujJm{uj) exp{iu!t)/[exp{f3uj) — 1]. For our present 
purpose of illustrating protection against SDE, it suf- 
fices to assume ohmic spectral densities with the same 
cutoff frequency ojc, that is, Jm{'^) = ?7mW exp(— w/wc), 
where rj„i, for m = 1,2,3, are dimensionless constants 
giving the respective strengths of dephasing, bit flip- 
ping, and dissipation. Calculating the continuum ver- 
sions of 1Cm{t) and Cm{t), using the ohmic spectral den- 
sities, gives K.m{t) = Tjm'jji/ {1 — i^ctf' and Cm{t) = 
(?7™//32)^'(i) (1 l/(/3wc) - itjli), where is the first 
polygamma function. By substituting these results into 
Eqs. ([H]), ([15|), ([H]), ([17|), and ([18|), we obtain the corre- 
lations that appear in Eq. (jlip . where the rotation ma- 
trix elements are obtained from Eqs. ([3|), and ([HJ. 
Once we have the coefficients Dp^q{t), for p,q = 1,2,3, 
then we can solve Eq. ([T2|) numerically. 

Here we impose the extreme situation where the disen- 
tanglement occurs faster than the time scale defined by 
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FIG. 2: Fidelity as a function of time for Bell states. The 
solid, dashed, dotted, and dot-dashed lines give the fidelities 
when the control fields are turned on, using N = 2,3, 5, 8, 
respectively. The insets (a) and (b) show, respectively, the 
fidelity and concurrence, both calculated att — t, as functions 
oiN. 

the inverse of the cutoff frequency, 1/ujc, and the thermal 
correlation time, tb = /S/tt. Hence, we take t = 27r/ti>c 
and, therefore, tc = t/N, as discussed below Eq. ([5]). 
For concreteness, in the numerical calculations we choose 
Tir^ — 2, Uz = 1, T — 10~^°s, a temperature of T = O.IK, 
771 = 1/16, 772 — 1/64, and 773 = 1/256. As initial pure 
states, we consider two classes: 

|<i>(0)) =cos0 ITT) + sine lii), (19) 
|*(0))=cos0|Ti)+sin0|iT), (20) 

for 9 £ [0, 27r), where we adopt the usual spin notation. 

To measure entanglement we use the concurrence 11], 
defined as the maximum between zero and A(i), with 

A{t) = Ai - A2 - A3 - A4, (21) 

where Ai > A2 > A3 > A4 are the square roots of the 
eigenvalues of the matrix p{t)ai^ya2,yP* {t)ai^ya2,y, where 
p*{t) is the complex conjugation of p{t), the reduced den- 
sity matrix of the two qubits in the Schrodinger picture. 
Thus, if A{t) gets less than or equal to zero, there is 
no entanglement and the state is separable. Our aim 
is to use external fields to keep A{t) fixed. Figure [T] 
shows A{t) for the evolution of different initial entan- 
gled states. In panel (a), A{t) is shown for initial states 
with A(0) w 0.70711, when the control fields are turned 
off (indicated with = 0) and on (A^ = 1,3,8). The 



solid and dotted lines correspond to the initial states of 
Eq. (Ull), with 6* = 7r/8 and 6* = Stt/S, respectively. The 
double-dot-dashed line corresponds to the initial states 
of Eq. (201), with 6* = vr/S and 6* = Stt/S giving the same 
A{t). In panel (b) A{t) is given for Bell initial states, for 
which A(0) ~ 1. The initial states of Eq. (fTO|) evolve 
resulting in the same A{t) for 9 = ±7r/4. Also, the initial 
states of Eq. pp]) evolve resulting in the same A{t) for 
6 = ±7r/4. In this figure panel, we show the results for 
Eq. p^ . since they differ from those of Eq. (^0]) only 
in the fourth decimal place. The solid line represents the 
result for the control field turned off, while the dotted, 
double-dot-dashed, and dashed lines show the results for 
control fields with A^ = 1, 3, 8, respectively. 

From the general theory of dynamical decoupling 
the protection gets better as N gets larger. Figure [2] 
shows the fidelity, F{t) = TT[pj{t)p{0)], as a function of 
time for Bell initial states. All four Bell initial states 
result in the same fidelity function up to the fourth dec- 
imal place. The solid, dashed, dotted, and dot-dashed 
lines give the fidelities when the control fields are turned 
on, using A^ = 2,3,5,8, respectively. The insets (a) and 
(b) show, respectively, the fidelity and concurrence, both 
calculated aXt — t, as functions of A^. 



V. CONCLUSION 

In summary, we have shown that it is possible to pro- 
tect a two-qubit entangled state from disentanglement, 
using a simple combination of a static field along the 
X axis and a rotating field in the yz plane. We have 
tested the method under a very unfavorable circum- 
stance, where dephasing, bit fiipping, and dissipation, at 
a finite temperature, are so effective as to disentangle an 
unprotected state within a time interval shorter than the 
characteristic reservoir correlation time, 2'it/uJc, and the 
thermal correlation time, tb = P/tt. Even so, the concur- 
rence can be preserved at high fidelity, as shown in Figs. 
[T] and [51 The present result also suggests that it might 
be possible to protect against disentanglement during the 
execution of an entangling quantum operation. 
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